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Redshift-Space Distortions with the Halo Occupation 
Distribution II: Analytic Model 
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ABSTRACT 



We present an analytic model for the galaxy two-point correlation function in 
redshift space. The cosmological parameters of the model are the matter density flm, 
power spectrum normalization erg i a-nd velocity bias of galaxies a^ , circumventing the 
linear theory distortion parameter (3 and eliminating nuisance parameters for non- 
linearities. The model is constructed within the framework of the Halo Occupation 
Distribution (HOD), which quantifies galaxy bias on linear and non-linear scales. We 
model one-halo pairwise velocities by assuming that satellite galaxy velocities follow a 
Gaussian distribution with dispersion proportional to the virial dispersion of the host 
halo. Two-halo velocity statistics are a combination of virial motions and host halo 
motions. The velocity distribution function (DF) of halo pairs is a complex function 
with skewness and kurtosis that vary substantially with scale. Using a series of col- 
lisionless N-body simulations, we demonstrate that the shape of the velocity DF is 
determined primarily by the distribution of local densities around a halo pair, and at 
fixed density the velocity DF is close to Gaussian and nearly independent of halo mass. 
We calibrate a model for the conditional probability function of densities around halo 
pairs on these simulations. With this model, the full shape of the halo velocity DF 
can be accurately calculated as a function of halo mass, radial separation, angle, and 
cosmology. The HOD approach to redshift-space distortions utilizes clustering data 
from linear to non-linear scales to break the standard degeneracies inherent in pre- 
vious models of redshift-space clustering. The parameters of the occupation function 
are well constrained by real-space clustering alone, separating constraints on bias and 
cosmology. We demonstrate the ability of the model to separately constrain Qm , erg , 
and a^ in models that are constructed to have the same value of /3 at large scales as 
well as the same finger of god distortions at small scales. 
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1 INTRODUCTION 

The growth of structure through gravity creates pecuhar 
velocities with respect to the smooth Hubble flow, mak- 
ing the line of sight a preferred direction when redshift 
is used as a measure of distance. The systematic differ- 
ences between the distribution of structure in real and red- 
shift space contain information about both the cosmology 
of the univer s e and the bias of galaxies observed ( see, e.g., 
iPeebledlToT^ : ISargent fc TurneiJllQTTl: iKaiseJIlQS'iD . In Pa- 
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per I JTinker et al.l2006cJl . we used numerical simulations to 
develop a blueprint for determining cosmological parameters 
from redshift-space distortions. In this paper, we use those 
simulations to calibrate an analytic model for redshift-space 
galaxy clustering that has the necessary accuracy to inter- 
pret high-precision measurements from large-scale galaxy 
redshift surveys, l ike the Two-Degree Field Galaxy Redshift 
Survey (2dFGRS: ^olless et al|l200lll and the Sloan Digital 
Sky Survey (SDSS: lYork et aLlboOOH . 



The effect of redshift-space distortions is most apparent 
at small scales, where the internal motions of galaxies in 
groups and clusters spread out close pairs along the line of 
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sight, creating the so-called "fingers-of-god" (FOGs). This 
picture of clustering in redshift space can be expressed as a 
convolution of the real-space two-point correlation function 
^Rij) with the probability distribution function (PDF) of 
galaxy pairwise velocities P{v^), i.e., 



l + C('-<T,r-„) = 



[l + U{r)]P{v.)dv., 



(1) 



where (,R{r) is the real-space correlation function, r^ is 
the projected separation, r^ is the line-of-sight separation, 
r^ = r^ + z^ , and v^ = H{r^ — z), where H is the Hubble 
constant. E quation ([T| is often referred to as the "stream- 
ing model" JPeebleJlOSCf) . Although this approach has been 
used prim arily to model observations of s mall-s cale distor- 
tions (e.g., |PeeblejlT979t iDavis fc Peeblejll983t iBean et all 
Il983^. equat ion Q is vali d in t he linear regime as well 
JFisherlllQQST) . Scoccimarro (l2004f) demonstrated that, pro- 
vided P{vz) is correct, the streaming "model" is a valid de- 
scription of the relation between the real- and redshift-space 
correlation functions on all scales. 

As presented, equation Q contains no cosmology. The 
aforementioned investigations utilized the streaming model 
to estimate the velocity dispersion of galaxy pairs, which 
was used to estimate the matter density parameter Qm 
through a "cosmic virial theorem". More recent studies of 
redshift-space distortions have utilized a modified linear the- 
ory model of anisotropics that, when applied to the real- 
space galaxy power spectrum Pfl(fc), takes the form 



Pz{k, fi) = Pfl(fc)(l + /?m')'(1 + k''alfi''/2y\ 



(2) 



where f3 — Q^ /bg, bg is the linear bias parameter, and /i 
is the cosine of the angle between the wavevector k and the 
line of sight. The te rm (1-1- (3fi^)'^, derived from linear the- 
ory bv lKaiseil (|l983), models the coherent flow of matter out 
of underdense regions and into overdense regions. The last 
term on the right hand side represents an exponential dis- 
tribution of random, uncorrelated peculiar velocities, which 
dominates Pz on small scales and is meant to encapsulate 
the FOG effect described above. Equation ^, commonly re- 
ferred to as the "dispersion model" , has two free parameters , 
(3 and the galaxy velocity dispersion ah . IScoccimarrd 1120041) 
points out several deficiencies in this model, both in the 
inability of linear theory to properly describe anisotropies 
even in the large-scale limit, and in the oversimplification 
of using a single parameter at, which has no clear physi- 
cal definition, to model small-scale velocities. Consequently, 
equation ^ introdu ces a 10 — 15% syste matic error in the 
determination of /3 JHatton fc Colelll999l : Paper I), a level 
of error significant compared with the precision achievable 
with SDSS and the 2dFGRS. 

The goal of this paper is to create an analytic model 
for the redshift-space correlation function by combining 
the streaming model of equation Q wi th the Ha l o Oc- 
cupation Distribu tion (HO D; see, e.g., Jing et al. 1998; 
Ma & Frv 2000; Peacock fc Smith 2000; Seliak 2000; Benson 



20011; IScoccimarro et alj 1200 ll: iBerUnd fc Weinberg! 120021: 
Coorav fc Shethll2002ll . The HOD quantifies bias on both 
linear and non-linear scales for a given galaxy sample by 
specifying the probability P{N\M) that a halo of mass M 



contains A'^ galaxies of a given type, together with any spa- 
tial and velocity biases between galaxies and dark mat- 
ter within individual halos. The HOD has been utilized 
to model the real-space clustering of galaxie s in the SDSS 
JZehayi et alJl2004l2005l:lTinker et alil2005ll and t he 2dF- 



GRS llTinker et alj|2006bl; see also lvang et alJl2003ll . In this 
paper we extend the HOD model from real space to redshift 
space by providing a model for P{vz) which is physically 
motivated and empirically calibrated on numerical simula- 
tions. 

Several recent papers have presented calculations of 
redshift-space distorti ons using ha l o models of dark matter 
and galaxy clustering JSeliakll200ltlWhitell200il:lKang et all 
120021; ICoorav1l2004l:1skibba et alJl2005l ; ISlosar et alJ^OOfiTl . 
providing insight into the role of non-linear dynamics and 
non-linear bias in shaping clustering and anisotropy. How- 
ever, these studies rely on the same linear the ory compo- 
nent of equation (|5| for large-scale anisotropies. iKang et alJ 
(2002) show that their model only reproduces the dark mat- 
ter Pz{k,y?) from N-body simulations after introducing a 
(Jk parameter for the halos, even after the viri al motions of 
partic les within halos were taken into account. ISkibba et alJ 
i2005l) demonstrate the difficulty is modeling redshift-space 
galaxy clustering in the transition region between quasi- 
linear and fully non-linear regimes using a linear theory de- 
scription of halo velocities. Other recent papers have used 
the halo appro ach to model galaxy and dark matter ve- 
locity sta tistics fS heth et ai]|2(301al ; ISheth fc Diaferidl200ll : 
ISheth et al. 2001b). While the model outlined in these pa- 
pers is derived from first principles, in contrast to the cali- 
brated model presented here, it is still based on linear theory, 
which does not provide the required accuracy for a robust 
implementation of equation Q. The purpose of our model 
is less as a first- principles derivation of P{vz) than as a tool 
to extract information from forthcoming observational data. 
In this context, the accuracy of the model is the paramount 
concern. In the course of developing the model, we will also 
gain new insight into the physics that determines P{vz), es- 
pecially the role of environment in producing a non-Gaussian 
velocity distribution. 

An accurate model for ^(rcr,T-^) with the HOD must 
properly incorporate halo motions. A proper model for halo 
pairwise velocities Ph{v) must correctly describe the distri- 
bution function for an arbitrary pair of halo masses, at any 
angle with respect to the line of sight, and as a function 
of separation. In the large-scale limit, linear theory is ade- 
quat e for describing the mean infall velocities of halos (see, 
e.g., Juszkicw icz et alJll999l ; ISheth et al. 2001a). However, 
the applicability of linear theory is problematic at scales 
where the observational data are robust. At all scales, linear 
theory docs not accurately predict the pairwise dispersion 
(Scoccimarro 2004). Higher order moments also play an im- 
portant role in Ph{v). N-body results have shown that the 
radial velocity PDF of dark matter halos exhibits significant 
skewness and kurtosis JZurek et al.lll994l ; Ijuszkiewicz et alJ 
ll998r .^ The skewness arises from the infall of matter into 
overdense regions (.Juszkiewicz ct al.. ■199&'l . The kurtosis, 

^ In this paper we use the convention of radial being the direction 
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manifesting as exponential wings in both the radial and tan- 
gential velocitie s, is du o to local non-linear effects for each 
halo in the pair. IScoccimarro ( 2004") concludes that a Gaus- 
sian is never a good description of velocities, even at the 
largest scales. IScoccimarrd i2004l focuses on velocity statis- 
tics of dark matter, but local non-linear effects apply to halos 

J II riK 

as well l|Kang et al.M20n2ri . and Ph{v) from simulations are 
non-Gaussian at all scales. It is not sufficient for Ph(v) to 
describe the first two moments of the velocity distribution. 
To accurately model C(''CT,f^), Ph(v) must reasonably de- 
scri be higher order mom ents of the distribution as well (see, 
e.g.. lFisher et al.l ll994V 

As stated in Paper I, our method for analyzing redshift- 
space distortions is to first use measurements of the pro- 
jected correlation function Wp{rp) to determine the param- 
eters of the HOD for a given cosmology. If HOD parameters 
cannot be found that allow the cosmological model to re- 
produce the observed Wp{rp) then that model is ruled out. 
Once the HOD has been determined, the redshift-space clus- 
tering is investigated by the analytic model presented here or 
the N-body approach of Paper I. The cosmological param- 
eters that most directly inffuence redshift-space clustering 
are the matter density parameter n^ , the amplitude of the 
linear matter power spectrum, defined here by erg, the rms 
linear-theory mass ffuctuation in 8 h~^Mpc spheres (where 
h = Ho/wo kms~^ Mpc~^), and the velocity bias of the 
galaxy sample, which we parameterize by a^, the ratio be- 
tween the satellite galaxy velocity dispersion and the virial 
velocity dispersion of the dark matter halo. One can con- 
sider models in which the velocities of central galaxies are 
biased with respect to mean motion of the host halo. In Pa- 
per I we considered the shape of the linear matter power 
spectrum Piin(fc) to be determined by measurements of the 
cosmic microwave background (CMB) anisotropies and the 
large-scale galaxy power spectrum (e.g., | Percival et alj2002l : 



ISpergel et alj|2003l: iTegmark et al.ll2004all . An advantage of 
the analytic model for ^(ro-, )",r) is that it allows for marginal- 
ization over Piin(fc) without the cumbersome process of run- 
ning multiple N-body simulations. The redshift-space ob- 
servables explored here and in Paper I are largely indepen- 
dent of Piin(fc), thus the constraints on HOD parameters 
come primarily from Wp(rp) and to a small degree the as- 
sumed Piin(fc), and the constraints on our cosmological pa- 
rameter space of Qrn, OS, and a„ come from the redshift- 
space data. Any degeneracies between cosmology and HOD 
parameters can be identified and marginalized over as well. 
In §2 we outline the analytic model, presenting analytic 
expressions for ^{rcr,rTr) by combining the HOD model for 
real-space clustering with the streaming model of equation 
Q. In §3 we detail the model for halo pairwise velocities, 
the key ingredient in the two-halo term, and calibrate it on 
the N-body simulations. §4 tests the model against the mock 
galaxy samples of Paper I, demonstrating that the model can 
recover the correct cosmological parameters from a wide va- 



connecting the halo pair, tangential being in a direction orthog- 
onal to the radial, and line-of-sight to be the direction from the 
observer. Velocities in these directions will be referred to as Vr, 
vt, and Vz, respectively. 
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Figure 1. A slice through one simulation showing the locations 
of halos of different masses. Large open circles are halos of mass 
M > 3 X 10^''. Small filled circles are halos of mass ~ 2 X 10^^ 
h~ Mq. The inset box plots the probability distribution function 
(PDF) of radial velocities for the lower-mass halo pairs in the 
slice. The shaded regions of the histogram highlights the skewness 
toward high- velocity pairs. All halo pairs within the shaded region 
of the histogram are connected in the slice plot. 



riety of cosmologies which produce significant degeneracies 
in redshift-space clustering. In §5 we summarize and discuss 
prospects for applying the model to current data sets. 



2 MODELING THE CORRELATION 

FUNCTION IN REAL AND REDSHIFT 
SPACE 

The real-space galaxy correlation function 5ij(r) can be cal- 
culated analytically with the HOD for a given galaxy sample 
and cosmology. In the HOD framework, galaxy clustering is 
separated into two distinct parts; pairs of galaxies within 
a single halo and those from two separate halos. The total 
correlatio n function is the sum of these two contributions 
fsee. e_.g..lBerlind fc Weinberdl2002l : ICoorav fc Shethll2002l : 
lZheng'2004"; and lTinker et al.l2005l for details on calculating 
the real-space correlation function). 

The real-space one-halo term is written as 



1 + Cih(r) 



2nr^n^ 



dM 



dn {N{N-1))m 
dM 2 



'2i?vir(M) 



F'(^—\ 



(3) 



where fig is the mean number of density of galaxies, dn/dM 
is the halo mass funct ion (here we use the formula presented 
in l Jenkins et al.l2001f) . and {N{N—1))m is the average num- 
ber of galaxy pairs in a halo of mass M. The function F' (x) is 
the fraction of galaxy pairs at the scaled radius x = r/2i?vir. 
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In practice, F'{x) is different for gafaxy pairs wfiicfi invoive 
the centrai galaxy and tliose between two satellite galaxies, 
and the calculation is separated into these two terms. 

In real space, one-halo pairs dominate £,R{r) at small 
scales. In redshift space, these pairs have large relative mo- 
tions due to the virial dispersion of the halos they oc- 
cupy, spreading these pairs out along the line of sight. 
We model the velocity distribution within each halo as 
an isotropic, isothermal Gaussian distribution, an approx- 
imation supported by numerical hydrodynamic simulations 
( Faltenbachcr et al. 2004) and observational analysis of rich 
SDSS clusters (McKay et al., in preparation). The satellite 
galaxy velocity dispersion in a halo of mass M is propor- 
tional to the virial dispersion. 






(4) 



^ ccn 

C 2 I 2 



We define i?vir to be the radius at which the mean inte- 
rior density of the halo is 200 times the background density, 
R^i^ = {4:Mh/3np200y/^. The constant of proportionality 
Uv is the galaxy velocity bias. For satellite pairs, the dis- 
persion is ass = \/2crsat. For central galaxies, a natural as- 
sumption is that these galaxies are at rest with respect to 
the center of mass of the halo, or a^c = 0. In Paper I we 
tested this assumption, demonstrating that Uvc = 0.2 has 
a negligible effect on redshift-space observables. We define 

— Q:„cO"vir, SO the disperslou of central-satellite pairs is 
' avir. Deviations from isotropy have lit- 
tle effect on ^(ro-,r^) for any physically reasonable value of 
the anisotropy, while non-isothermality resulting from spa- 
tial bias or unrelaxed systems can be modeled as an "effec- 
tive" velocity bias parameter (Paper I). 

In equation ^, the total number of one-halo pairs in- 
volves an integral over the halo mass function. At a given 
separation, the velocity PDF for one-halo pairs is a super- 
position of Gaussians weighted by the relative number of 
galaxy pairs from each halo of mass M. This superposi- 
tion of Gaussians of varying dispersion makes the overall 
pair wise galaxy v elocity dispersion approximately exponen- 
tial_JShBth[^96f) , in a ccord with observat ional inferences 
JPavis fc Peeb les 1983; iHawkins et al.ll20o3) . 

The combination of equations Q and @ for central- 
satellite pairs is 



^[fir^^r^) 



^^1^1 Jo 



dM 



dn {Nss,t) M (Nccn) M 



dM 



F' 



, / Vr'i + z^ 



X exp 



2itvir 

(ry ^ zf 
2ols 



2-Rvir 

1 



27rCTc 
dz 



(t1 + 22) ' 



(5) 



where {Ncon)M and {A'aat)M are the mean number of cen- 
tral and satellite galaxies for mass M, respectively. _F^j,(x) 
is proportional to the halo density profile, which we assume 
to follow the form of Navarro et al. (1997, hereafter NFW). 
For satellite-satellite pairs, the one-halo term is 



^[f{r„,r^) 



dM 



dn {N.^tiNs, 



dM 
Vr'i + z' 



X exp 



(r2+22)' 



(6) 



where {Nsat{Nsat — 1)) is the second moment of the satellite 
occupation functio n, and Fss(x) is the fraction of satellite- 
satellite pairs (see ISheth et al.l l2001al for a derivation of 
Fcs(x) and Fss{x) for the NFW profile). The full one- halo 
correlation function is ^ih{rcr,r,T) ~ ^J" -I- Cih ■ 

The two halo term is most straightforward to calculate 
through a direct implementation of the equation Q , 



1 + C2h(r<„r^)= / [l + ^i,{r)]P2bMr,,p)dv, 



(7) 



where C2h('') i^ the two-halo contribution to the real-space 
correlation function, P2h(yz\r,(t>) is the line-of-sight velocity 
PDF of galaxy pairs from two distinct halos, and cos = 
Vcr/r. This PDF is a convolution of the line of sight PDF 
of halo pairwise velocities with the Gaussian distribution 
for galaxy velocities within each halo. P2h{vz\r, <j)) is a pair- 
weighted average integrated over all possible combinations 
of halos; 



P2h(t'.|r-,0)=2(n; 



Mii, 



rfM,_(iV)M, 



dn 



dM2-^^{N}M,Ps+i,{v4r,<l,,Mi,M2) , (8) 



where Pg+h is the velocity PDF for galaxy pairs from ha- 
los of masses Mi and A/2. Pg+h includes both the internal 
motions of the galaxies within each halo and the relative 
motions of the halo centers of mass. The limits on each inte- 
gral in equation JSJ are determined by halo exclusion; halo 
pairs cannot be closer than the sum of their virial radii (i.e., 
Rvii,i + Rvii,2 ^ r). At large separation, this effect is neg- 
ligible because the largest collapsed objects have radii ~ 2 
h~^Mpc, and the limits of equation JSJ approach infinity 
when r is large. At small scales, halo exclusion strongly in- 
ffuences the number of pairs. In equation ©, Miim.i is the 
maximum halo mass such that -Rvir(Miini,i) — r— -Rvir(A^miii) 
and Miiin,2 is related to Mi by -Rvir(Miini,2) — r — Rvii{Mi), 
where Mmin is the minimum mass halo that can host a 
galaxy. The restricted number density, n'g , is the total num- 
ber of two-halo galaxy pairs at separation r, also calculated 
using halo exclusion by 



dMi 



dn 
dM 



{N)mi 



A-fii„ 



dMa 



dn 



{N)m,. (9) 



We denote the velocity PDF of halo pairs as 
Ph(vz\r,(f),M\,M2)- This represents the center-of-mass 
motions only. Once the relative motions of the dark matter 
halos are determined, Pg+h is calculated by convolving Ph 
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Figure 2. Velocity PDFs for halo pairs at different local densities. The top four panels plot PDFs as a function of local environment 
for low- mass halos (solid histograms) and high-mass halos (filled circles connected with dotted lines). The local density is calculated by 
a top-hat smoothing kernel centered on the midpoint between the halo pair, with radius equal to the halo separation. The bottom panel 
shows the unconditional velocity PDF for the same halos. 



with the Gaussian describing the internal motions of the 
galaxies within the halos. The model for Ph{vz\r, 4>, Mi, M2) 
is important both for capturing the behavior of the redshift- 
space correlation function in the transition region between 
the one- and two-halo terms and for correctly modeling dis- 
tortions between quasi-linear to linear scales. We describe 
it in detail in the next section. 



3 THE HALO VELOCITY MODEL 

3.1 Conditional matter PDF 

In §1 we listed the requirements an accurate model of Phiv) 
must satisfy. To meet these requirements, we use the N- 
body simulations from Paper I to explore halo velocity 
statistics. Our simulation set consists of five realizations 
of an inflationary ACDM cosmology with Q.rn =0.1 and 



(Tg = 0.95 at 2: = 0, using the public N-body code gad- 
get iSpringel. Yoshida. fc Whit3l200J) . We utilize outputs 
at higher redshifts to represent models at different values 
of (jg and Q,m- Each simulation is 360'' particles in a cubic 
volume 253 /i~^Mpc on a side, yielding a mass resolution 
of 9.64^" X Q,m h~^yiQ. The simulations have a power spec- 
trum shape parameter F = 0.2 in the parameterization of 
lEfstathiou et all (|l992), and a spectral index n^ = 1. See 
Paper I for further details regarding the simulations. Our 
fiducial cosmology is Q.m = 0.3, (Tg = 0.8, corresponding to 
z = 0.55. Unless otherwise stated, N-body results will be 
using this output. 

Figure elucidates the complexity of 
Ph{v\r,4>,Mi,M2). The inset box shows the PDF of 
radial velocities of halo pairs of mass Mi = M2 ~ 2 x 10^^ 
/i~^M0 at a separation of r ~ 10 /i~^Mpc. Negative radial 
velocities are defined as the halos going toward one another. 
The PDF is not well described by either a Gaussian or an 
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Figure 3. Comparison between the conditional matter PDFs measured from the simulations {solid histograms) and the model of 
Equation 1111 (dotted lines). Top row: Pm{&\i',Mi,M2) for fixed mass {10-*^^ h~^y[Q) and erg = 0.8, but changing radius. Middle row: 
^m(5|r, Mi,M2) for fixed radius (10 /i~^Mpc) and eg = 0.8, but increasing halo mass. Bottom row: Pm{5|r, Mi, M2) for fixed mass 
(lO'^^-^ h-'^Mo) and radius (10 /i-^Mpc), but changing erg. 



exponential. The shaded regions of the PDF highlights the 
skewed tail of halo pairs rapidly approaching each other. 
The outer panel plots the positions of these halos in a 40 
/i~^Mpc slice through one realization. Small, filled circles 
indicate the 2 x 10^^ /i~^M0 halos of the inset histogram. 
The open circles in the plot represent high mass halos, 
M > 3 X 10^^ /i"^Mq. AU pairs inside the shaded region 
of the histogram have been connected with a line. These 
high-velocity pairs are almost exclusively in regions that 
include a high mass halo. Because high mass halos prefer- 
entially occupy dense regions, we recover the well-known 
result that the non-linear halo velocity field is coupled to 
the non- linear density field. However, Figure Q leads us to 
investigate the velocity PDF as a function of environment. 
In Figure |5| the lower panel plots the same PDF from 
Figure Q for 10^^ fe-^Mg halos at r ~ 10 /i"^Mpc. This 
panel also plots the PDF for 3 x 10^^ H'^Mq at the same 
separation. The width of these mass bins is a factor of 



two. The PDF for high-mass halos is noisier due to lower 
statistics, but it is apparent that the PDF for high-mass 
halos peaks at u ~ —150 kms~^, while the mode of the 
low-mass PDF is at zero. It is unlikely to find high-mass 
halos moving away from one another at this scale, while 
that probability is nearly 40% for low-mass halos. 

The upper panels of Figure |21 show the PDFs for the 
same halos, now binned by local overdensity, 5 = p/p—1. We 
calculate 5 inside spheres of radius equal to the separation 
of the halo pairs, centered on the midpoint between the two 
halos. Setting the radius of the smoothing kernel equal to the 
halo separation allows for cleaner determination of the pair's 
environment. The four top panels span a range in mean local 
density of (5) ~ to 3. From these data the origin of the 
high-velocity tail is clear; the negative skewness arises from 
halo pairs which lie in dense environments. At 5 = 0, the 
velocity distribution is narrow and peaked at u ~ kms~^. 
As S increases, both the mean velocity and the dispersion 
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Figure 4. The scale and mass dependence of parameters of the 
two-halo velocity model. Each line represents of halo pair with 
mass ratio M2/M1 = 4% with i = - 3 and Mi = lO^^ h-^MQ. 
Panel (a): cutoff density scale for the 2-halo conditional mat- 
ter PDF (Eq. 1131 ). Panel (b): Power-law index for the density- 
dependence of the pairwise velocity dispersion (Eq. I19I V There 
is only one line in this panel because there is no dependence on 
mass. Panel (c): Normalization of the density dependence of the 
pairwise radial velocity dispersion (Eq. 1211 '). Panel (d): Normal- 
ization of the density dependence of the pairwise tangential ve- 
locity dispersion (Eq. 1221 '!. 



increase as well. At fixed 5, the PDFs no longer show strong 
skewness or kurtosis; a Gaussian is a good approximation. 
Therefore the total distribution in the bottom panel can be 
modeled as a superposition of Gaussians weighted by the 
probability of finding a halo pair at a given 5. 

At a given density, Figure |5| also demonstrates that the 
PDFs for high- and low-mass halos are very similar. It is 
the different distributions of local environments which leads 
to the differences in the total PDFs in the lower panel. We 
therefore propose the ansatz 



Ph{vir,t]\r,Mi,M2) = 

PG(v[r,t]\r, S, Ml, M2) Pm{S\r, Mi , Ma) dS, (10) 

to model the radial and tangential pairwise motions of ha- 
los. In equation llUl Pq is a Gaussian, and P,n{S\r, Mi, M2) 
is the conditional matter density PDF in spheres of radius r 
centered on a halo pair of masses Af 1 and M2 . An ansatz of 
this type has been used to model one-point velocity statistics 
of ha los and galaxies ( Shetli fc Diaferio 2001 ; Hamana et aJj 
l2003l) . Ahhough the PDFs in the upper panels of Figured 
appear independent of halo mass, Pa varies moderately with 
niEiss at small scales because at small r the halos themselves 
contribute a large fraction of the total local mass. At scales 



significantly larger than the halo radii, Pq becomes inde- 
pendent of mass. 

The PDF for the unconditional matter distribution 
Pm{S\r) is well fit by a lognormal distribution function 
with dispersion that varies with smoo thing scale (e.g., 
[Colcs & Jones 1991; Kofina n et al.lll994) . The relationship 
between the halo mass function and t he large-scale den- 
sity field can also be quantified (e.g., jBonde t alJ Il99ll : 
ISheth fc Tormenll200d iPavlidou fc Field jl2005^ . The con- 
ditional mass function can be roughly approximated as 
n{M\S) ^ [1 + b{M)S]n{M) (fMo fc White 1996), where 
6(M) is the large-scale bias factor for halos of mass M (here 
we use the halo bias formula in Appendix A of lTinker et alJ 
I2OO5). The distribution of densities at which a halo of a 
given mass can be found is roughl y Pm{5\r,M) « [1 -|- 
fe(M) S] Pm{S\r) iHamana et alj2003i) . The term [l + b{M) S] 
effectively shifts the matter PDF to higher densities by an 
amount which depends on the bias of the halo. 

The density distribution around pairs of halos, de- 
noted as Pm(5|r, Afi, M2), is a murky theoretical prob- 
lem with no obvious straightforward approach. To model 
Pm((5|r, Ml, M2) we truncate the unconditional matter dis- 
tribution with an exponential cutoff at low densities, i.e., 



P™((5|r,Mi,M2) = Aexp 



po(r-,Mi,M2) 



(11) 



where A is a normalization constant to make the total prob- 
ability unity, p = 1 + 5, and po is a density cutofi' scale to be 
calibrated from the simulations. The one-point non-linear 
distribution function for the dark matter takes the form 



'\r) 



2ttoI 



exp 



:in(l + S)+ ai/2f 



2af 



(12) 



Om jr) is the mass vari- 

;e.g.,ir 



where o-\{r) — ln[l + <y'm{r)] and 

ance in top-ha t spheres of radius r (e.g.. lColes fc Jonesll99ll : 
iKofman et al.l | l994,'l . To calculat e amjr), we u se the non- 
linear matter power spectrum of ISmith et al.l i2003r) . The 
cutoff density scale po is 



po^pi [h{Mi) + h{M2)]+(—\ 

\roJ 



(13) 



At large scales, equation I13II is proportional to the sum of 
the halo bias factors. At small scales, where the masses of 
the halos themselves contribute significantly to local density, 
the mass contained within r cannot be below M1+M2. Thus 
the cutoff density should scale roughly as po oc (r/i?vir)~''. 
From our N-body simulations, we find pi = 1.41, ao = —2.2, 
and ro = 9.47?vir,i, where 7?vir,i is the virial radius of the 
larger halo of the halo pair, provide a good fit to the data. 
Figure 121 compares equation lllll to Pm{S\r,Mi, M2) mea- 
sured from the N-body simulations. The top row of pan- 
els plots log[P((5)] measured in spheres centered on halo 
pairs of separations 5, 10, and 20 /i~^Mpc, for halos of 
mass 10^^ h~^MQ. The N-body results, shown with the solid 
histograms, are lognormal in shape like the unconditional 
Pm{S), with a dispersion that narrows as the smoothing 
scale increases. Dotted lines show equation IIH . with po 
calculated from equation 11311 using the halo bias function 
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Figure 5. Comparison between the velocity statistics of N-body 
halos and the velocity model integrated over all halo masses. In 
all panels the solid line plots the N-body results and the open 
circles connected by dotted lines plots the model. The top four 
panels present the first four moments of the radial velocity PDFs 
as a function of scale: mean /i, dispersion Cr, skewness /^s/cr^, 
and kurtosis k. The bottom two panels show the dispersion at 
and kurtosis k of the tangential PDFs. For proper comparison, 
the integrals over velocity for the moments are binned in the same 
manner as the N-body results and are truncated at the high and 
low velocities of the N-body statistics. 



of iTinker et alJ (|200a). Equation IIH accurately tracks the 
change in Pm{S\r,Mi,M2) with smoothing scale. 

The middle row of panels in Figure compares the 
model to the N-body results for different halo masses at the 
same smoothing scale of 10 h~^Mpc. N-body halo pairs of 
mass ~ 2 X 10^^ h~^MQ, which have large-scale bias of ~ 0.8 
in our fiducial cosmology, reside on average in overdense re- 
gions, with (5) = 0.80. As the mass increases by factors of 
4 and 16 (with b ~ 1.1 and 1.5, respectively), (S) increases 
to 1.00 and 1.36. The mean density calculated from equa- 
tion mil l for each halo mass is 0.82, 1.02, and 1.32, in good 
agreement with the numerical results. 

The bottom panel in Figure |21 isolates the effect of 
changing ug on Pm(5|r, Mi, M2). The three panels present 
numerical results for A/ ~ 5 x 10^^ h~^MQ halos at r ~ 10 
h~^Mpc for as = 0.6, 0.8, and 0.95. In equation (EU, the 
dependence of the matter distribution on the power spec- 
trum is through the parameter arn{r), which is proportional 
to as. As as increases, the dispersion in the unconditional 
Pm{S\r) from equation (I12II increases. The same is true for 
Pm{5\r, Ml, M2); the N-body results broaden substantially 
from as = 0.6 to 0.95. The results from equation lllll . plot- 
ted once again with the dotted hues, model the changes of 
Pm(5|r, Ml, M2) with as accurately. With three adjustable 
parameters, the quantities pi, ao and ro of equation 113II . 
our model for the conditional matter PDF Pm(S\r, Mi, M2) 



reproduces the distinct effects of varying halo mass, smooth- 
ing scale, and power spectrum normalization. 



3.2 Parameters of the Gaussian 

For a halo pair at separation r and density 5, we approximate 
the radial velocity PDF as a distribution of the form 



P{vr\S,r,Mi,M2) = 



2'JTar 



exp 



-{v- 



Mr)^ 



2af 



(14) 



Both the mean fir and the dispersion ar of the radial velocity 
PDF depend on 5. The PDF for tangential velocities is also 
a Gaussian, but with dispersion at and zero mean. At large 
scales, or at negative density contrast, /i,, is well described 
by linear theory, 

f,U5,r) = -Hrnt'l- (15) 

At small scales, Hr is approxi mated by the n on-linear spher- 
ical collapse model (see, e.g.. |Peacocklll999r . In this model, 
the dependence of the velocity perturbation on density is 
expressed parametrically as 



^ ^ 9 (7 - sin 7) _ ^ 

2 (1 — cos7)'^ 

3 sin 7 (7 — sin 7) 

(1 — cos 7)2 

The mean radial velocity is 
fi,,{5,r)^Hrnlfu{S) exp[-(4.5/rp)'] 



(16) 
(17) 

(18) 



where the exponential term is added to better match the 
N-body results at the smallest scales but has little influ- 
ence on the overall behavior of /isc. Investigation of our N- 
body results shows that the spherical collapse model best 
describes the N-body simulations at r < 4 h~^Mpc, while 
linear theory is an excellent description of the results for 
r > 20 h~^Mpc. In the transition region, we express /ir as a 
weighted mean of equations 11151 and 11811 . which smoothly 
transitions between the two regimes. The weighting fac- 
tor, w^ increases linearly from to 1 in ln(r) in the range 
4 ^ r ^ 20 /i^^Mpc, with fir = Wfisc -I- (1 — w)fiiin. 

Recalling Figure |5| the skewness of the radial veloc- 
ity PDF arises from the combination of halo pairs in high- 
density regions, which have a large streaming velocity and 
high dispersion, with halo pairs in mean- to low-density re- 
gions, which have low mean streaming velocity (relative to 
Hubble flow) and small dispersion. In our N-body simula- 
tions, the skewness of the radial velocity distribution does 
not monotonically increase as halo separation decreases. As 
halo pairs become close, of order twice the sum of their 
virial radii, the skewness decreases and the PDF becomes 
more symmetric. Because fisc diverges rapidly as S becomes 
large, the skewness of the model will always increase with 
decreasing r. To compensate for this effect, we enforce the 
condition that at r = 4_Rvir,i, fJ-r is held constant for all 5 at 
the value of the most probable 6 for that separation, remov- 
ing the skewness entirely. In the full halo-|-galaxy PDF, this 
change from skewed to symmetric distribution functions is 
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Figure 6. Comparison between line-of-sight velocity PDFs measured from the simulations and from the analytic model. Solid histograms 
show the N-body results and dotted lines plot the model calculations. The sixteen panels create a grid in which halo mass increases from 
left to right and separation increases from bottom to top. Halo pairs are all of equal masses and at an angle of 45° with respect to the 
line of sight. 



smoothed out in scale through the change in -Rvir,i with halo 
mass. 

It is clear from Figure |5| t hat the velocity disp ersion 
also depends on the local density. iHamana et al.l ll2003r show 
that the one-point velocity dispersion of halos is described by 
a power law, i.e. a^ oc (1 -1-5)". We find that this is also true 
of the pairwise velocities but, that both the normalization 
and index of the power law are functions of scale. The power- 
law index varies weakly with radius, and is modeled by 



a{r) 



35/i-iMpc 



(19) 



In linear theory the dispersion of halo space velocities 
scales with cosmology as Q^^jfag. We parameterize the halo 



velocity dispersions in terms of the overdensity p200 at which 
the dispersion is 200 kms~^ in our fiducial cosmology; 



""('-,*) 



\0.3J \0.8J 



P2OO {r,t) 



kmf 



(20) 



where P200 (r,t) is a function of Mi , M2 , and r. We find that 
the following formulas provide a reasonable description of 
the radial and tangential velocity dispersions in the simula- 
tions: 



P200,r(Ml,M2,r) 
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Figure 7. (.{ra,r-„) for HODS in Table 1. The contours plot 
constant amplitude of 5(rCT,r-n-), separated by powers of 1.7. The 
outermost contour is 1.7~^ = 0.07. Thick gray lines show the 
N-body results and thin black lines are the model calculation. 



1/2 



5.0i?. 



p2Qo,t{Mi,M2,r) = 



+ 



11.57?, 



1/2 



+ 0.50, 



(21) 



7.2i?. 



1/2 



+ 



12.6i?if'o , 



+ 0.48, 



(22) 



where i?vir,o = Rvir,i+Rvii,2- The units of r and Rvir in equa- 
tions 1211 1 and 12211 are h~^Mpc. The mass dependence again 
enters with the virial radii of the halos, but the dependence 
on halo mass is weaker than for the density distribution in 
equation I13I I . Figure |1] shows the radial dependence of the 
four model parameters which must be calibrated to the sim- 
ulations, po, a, P200,r, and p200,r, from equations CHli GHli 
12111 , and 12211 , respectively. The different lines in each panel 
are for halo pairs Mi and M2, where Mi = 10^^ /I'^M©, 
and M2/M1 — 4\ and i = — 3. At large scales, the depen- 
dence on halo mass becomes negligible except when it enters 
in the bias factor used in the cutoff of the density PDF. At 
these scales, the value of p2oo approaches a constant for both 
tangential and radial velocities indicating that the velocity 
dispersion becomes independent of scale and mass. 

Using this model, we integrate over all halo pairs and 
calculate equation JHJ for halos alone, without any contri- 
bution from galaxies (i.e. {N)m ~ 1). For this special case, 
equation Q becomes 



Pn{v\r) = 2{n'^y 



A/ii„ 



dMi 



dn 
'dWi 



dMo 



dn 
■dAk 



Ph{y\r,Mi,M2) 



(23) 



where v can refer to either Vr or Vt, and n'^ is the num- 
ber density of halos allowed by halo exclusion, in analogy to 
equation I^J. Equation l|23^ allows a direct comparison to 
our N-body simulations. Figure |3 shows the first four mo- 
ments of the halo velocity PDF as a function of scale for 
the radial velocities, and the first two even moments for the 
tangential velocities, for our fiducial cosmology. The N-body 
results are shown with the solid line and the model is the 
open circles connected by the dotted line. For calculating the 
velocity moments with the model, we use the model value 
at the center of each bin that has a non-zero probability in 
the N-body distribution and set the more extended wings of 
P{v) to zero. The contribution from the high-n wings makes 
little difference for the comparison of the first two moments, 
but it noticeably affects the comparison of the skewness and 
kurtosis to the N-body results. 

The mean and dispersion of the radial velocities are 
accurately described by the model. For the skewness, de- 
noted as fj-s/a^ with /is — '^{vi — fi)"^ / N , the N-body results 
demonstrate the previous statement that the asymmetry is 
minimal at large separation, reaches a peak at intermediate 
scales near 4 /i~^Mpc, and turns over to become nearly zero 
at r < 1 h~^Mpc. The model describes this behavior well, 
reaching a peak at the same location and reducing to zero 
at r = 1 h~^Mpc. The kurtosis, defined as k = fXi/a'^ — 3 
with fi4 — "^{vi^ fJ.)'^/N, approaches zero at large scales for 
both the tangential and radial velocities, and monotonically 
increases with decreasing r to a value of k « 7 at r « 1 
/i^^Mpc (for reference, k = 2 for an exponential distribu- 
tion) . The kurtosis is moderately over-predicted at interme- 
diate scales in the model, but both the sign and slope of 
the trend with radius are correct. As we will show later in 
§5, this small discrepancy with the fourth moment does not 
adversely affect the accuracy of the (,{rcr,rTr) model. 



3.3 Ph{vz) as a function of angle 

For a halo pair at angle </!>, the line of sight relative velocity 
is Vz = Vr sin 0-1- «t cos (j), where cos = ra/r. Unfortunately, 
knowing the radial and tangential velocity PDFs indepen- 
dently does not allow one to compute P{vz) because the 
radial and tangential velocities are correlated; if the radial 
infall velocity is high, the tangential velocity is also likely 
to be high. However, this correlation is accounted for by the 
density dependence of the PDFs in our model, and at fixed 
5 the distributions P{vr) and P{vt) are independent. The 
line-of-sight distributions can then be computed from the 
appropriate convolution. 



Ph{vz\cj>,S,Mi,M2,r)^ 

Ph (vt I (t>, 5, Ml , M2 , r) Ph (vr 1 (p, 5, Mi,M2, r)dvt , 



(24) 
where Vr = {vz — «t cos(j>)/ siiKJ), to obtain 
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Ph{v,\(l>,S,Mi,M2,r) = 



1 / 2 2 I 2 • 2 i\~^l'^ 

at COS (p + a^ sm 01 



-(Vz — Sm(f) fj.r) 



X exp 



2(cr| cos^ (/) + (Tr sin 



(25) 

Inspection of equation 1251 produces the expected result 
that at = tan(^ = siiKJ) = 0, the expression reduces to 
a zero- mean Gaussian with dispersion at- For <j) — it/2, it 
reduces to equation I14I I. 

Rather than compare the moments of Ph{vz) from the 
model to the N-body simulations at a given 0, since they 
appear much as they do in Figure |K1 instead we compare 
the PDFs themselves for several masses and separations at 
an angle of 45°. The sixteen panels of Figure |S| represent a 
grid of models where each row is a different separation and 
each column is a different mass. In each column, the halo 
pairs are Mi = M2 = 4' Mc, with i = 1 - 4, and Mc = 10^^ 
h~^MQ. In each panel the solid histogram represents the 
N-body results and the dotted line plots the model calcu- 
lation. The agreement between the model and the N-body 
results is excellent wherever the N-body statistics are good. 
The model reproduces the numerical results over an order 
of magnitude variation in radius and almost two decades in 
halo mass. 



3.4 Combining halos and galaxies 

Because the halo PDF at fixed 5, r, and halo masses is a 
Gaussian, internal motions of the galaxies can be included 
by adding the galaxy dispersion to the halo velocity disper- 
sion in quadrature. As in the one-halo term, pairs involv- 
ing central and satellite galaxies must be treated separately 
once again. In the two-halo term, there are four distinct 
types of pairs with four velocity dispersions: central-central 
pairs, satellite-satellite pairs, and satellite galaxies from Mi 
paired with the central galaxy of M2 and vice versa. The 
four velocity dispersions are 



a^ = al + iaLcrvir(A^i) + ala^i,{M2)] 
al + [al,a^i^(M2) + Q^a^i,(Afi)], 



(26) 



-I 



where a^ = at cos^ (j) + a'^ sin^ (j). The complete PDF is a 
linear combination of these four Gaussian functions. 



Pg+hivz\(t>,r,Mi,M2,5) 



4 

E 

i=l 



27ro"i 



■ exp 



-(«. 



fir sm ( 



2a? 



(27) 



where each term is weighted by the relative number of pairs 
of each type, i.e., 



Wl = {iVccn)A/i {Nccn)M2/{N)M, (N) M^ 

W2 = (iVsat)A/l (A^sat)M2/(A^)Mi (N) M^ 

W3 = (iVccn)Mi (iVsat)M2/(A^)Mi (N) M2 

W4 = {Ns^.t)M^{Nccn}MJ{N)M^{N)M2■ 



(28) 



To obtain the full galaxy-|-halo PDF for a pair of halos 
Ml and M2, equation l|27t is weighted by P,„((S|r, Mi, A/2) 
and integrated over all densities, i.e.. 



Pg+h{vz\(t),r,Mi,M2) 



Pg+h{v4(t>, r, Mi,M2,5) Pm{S\r, Mi,M2)d5. (29) 



To complete the calculation of the two-halo velocity PDF, 
equation (I29II is integrated over all combinations of halo 
masses as in equation (|SJ , which gives P2h (^z |<^, r) (note that 
the subscript 2h refers to two-halo galaxy pairs, as in §2). 
This PDF is used in the two-halo streaming model, equation 
0. Once the total velocity PDF has been calculated and 
tabulated as a function of separation and angle, ^ih{ra,rTy) 
is calculated from equations and and ^2h(»'o-, J"^) from 
equation 0. 



4 RESULTS 

Although one-dimensional diagnostics are preferable when 
analyzing observational measurements of redshift-space dis- 
tortions, it is worthwhile to demonstrate the ability of 
the model to reproduce the overall features of the two- 
dimensional ^(ro-,r.^). Figure |7| plots ^(ro-,r.^) for both the 
N-body results and model for our central model, HODS in 
Table 1, plotted as contours of constant correlation in the 
To- — r^ plane. The thick gray contours represent the N-body 
data, averaged over the fifteen orthogonal projections of the 
five realizations, and the thin black lines plot the model re- 
sults. The contours are separated by powers of 1.7, with the 
outermost contour representing 1.7~^ = 0.07. The coher- 
ent infall of matter into overdense regions is present in the 
fiattening of the contours on large scales. In the absence of 
this infall, the contours would be concentric circles at these 
scales. At small scales the FOG effect is shown by the sharp 
elongation of contours at ro- < 2 h~^Mpc. The model ac- 
curately reproduces both the large- and small-scale effects 
seen in the N-body results. 

Although it is possible to compare model predictions 
directly to the observed two-dimensional correlation func- 
tion £,{rcr, r,r), this is not preferred due to the complexity of 
the error analysis. There is a strong covariance between data 
points, making it difficult to accurately estimate the covari- 
ance matrix of the measurements due to the large number of 
data points. Using the two-dimensional data set also makes 
it difficult to isolate constraints on the cosmological param- 
eters from constraints on the HOD parameters. We proceed 
instead by extracting one-dimensional observables from the 
two-dimensional data, observables that are more easily in- 
terpreted in terms of the parameters Qm, ag, and a„. 



4.1 Redshift-space observables 

In Paper I we described three observables for the correlation 
function. The first two, based on multipoles of ^(rCT,r^), 
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Table 1. Properties of the Simulations and HOD Parameters 



Model 


eg 


Af,„i„ [h-iMo] 


Ml [/i-IMq] 


OisaX 


!^m(/3 = 


.46) 


OLv 


HODl 


0.95 


1.12x1012 


2.81x1013 


0.934 


0.24 




1.109 


HOD2 


0.90 


1.12x10" 


2.75x1013 


0.959 


0.26 




1.076 


HODS 


0.80 


1.09x10^2 


2.51x1013 


1.005 


0.30 




1.000 


HOD4 


0.70 


1.08x10^2 


2.24x1013 


1.109 


0.36 




0.906 


HODS 


0.60 


1.06x1012 


1.93x1013 


1.199 


0.47 




0.797 



Note— All masses have been scaled to f2„ 



: 0.3. When we scale an HOD to a different value of ^m, the values of Mmin and Mi scale in 



proportion to f2m. Values in the sixth and seventh columns are the parameters of the HOD models used in §4.3, which have /3 = 0.46 

and the same finger-of-god strength for all five models. 



emphasized large-scale distortions because they have pre- 
dictions in linear theory. The multipoles of the correlation 
function are expressed as 



ii(r) = 



21 



i{r„,r^)Pi{p)d^i, 



(30) 



r~ = r^ + r^, ^ = r-^jr, and Pi (pi) is the Legendre 
polynomial of order I. The ratio of the monopole to the real- 
space correlation function, 



ioiR{r) = 



?o(0 



(31) 



is predicted in linear theory to be a constant value larger 
than one, but non-linearities in the velocity field suppress 
small-scale redshift-space clustering, driving this ratio to 
zero as r becomes small. In Paper I, we restricted our use of 
this diagnostic to the asymptotic value at r ^ 10 /i~iMpc, 
but small scales contain additional information, and in this 
analysis we will make use of the full range of r. 

The quadrupole moment of the correlation function is 
defined as 

6(r) 



Q«W 



^o(r-)-Co(r) 



where ^o('") is the volume averaged monopole, 



Ur) 



^o{s)s'^ds. 



(32) 



(33) 



In the numerical results of Paper I, this diagnostic also has 
an asymptotic value at large scales that is roughly consistent 
with its linear theory prediction, but due to the higher-order 
multipole, equation l|32|l is more sensitive to the galaxy ve- 
locity dispersion, and non-linearities affect this diagnostic 
at larger r than ^q/h- For cosmologies that produce strong 
redshift distortions, e.g. /3 > 0.6, Q^ does not reach a flat 
plateau even at 40 h~^Mpc, putting the linear theory regime 
outside the range of high precision measurements. 

To quantify small-scale distortions, Paper I focused on 
the behavior of ^{r-n) at flxed, small Vcr- We define the quan- 
tity r^/2 as the value of r^ at which the correlation function 
falls by a factor of two from its value at r^ = 0, or 



(34) 



m 2- 

The quantity r^/2 can be thought of as a characteristic 
length scale for FOGs. At r ^ 0.1 h~^Mpc, r^ji depends 



on the parameter combination fima^ and is nearly indepen- 
dent of (Tg- At intermediate scales, r ~ 1 /i~iMpc, rj/2 is 
highly sensitive to ag because the number of one-halo pairs 
at this scale depends on the shape of the halo mass func- 
tion at large masses. For the approach of Paper I, measuring 
r5/2 at Ta =0.1 hr Mpc is beneficial because it provides a 
simple relation to constrain QrnCe^- 

The form of the HOD used in Paper I is simple in its 
parameterization; halos above a minimum mass Mmin have a 
central galaxy located at the center of mass of the host halo, 
while satellite galaxies have a mean occupation number that 
scales as a power law with the host halo mass, i.e. (A'sat)A/ = 



(M/Mi 



Therefore the mean occupation function is 



{N)m = {iVccn)A/ + (iV., 






(35) 



for M ^ A/niin, and {N)m = below Afmin. The dispersion 
of satellite galaxies about {A'sat)M is a ssumed to be Poisson , 
motivated by both numerical results IIKravtsov et alJ jgOO^ 
[Zheng et al.l 12005^ and observational studies JLin et alJ 
120041) . The mean occupation of central galaxies is a step 
function with no scatter. This parameterization has only 
two free parameters^ but still has enou gh freedom to mode l 
the observed lOp(rp) of SDSS galaxies iZehavi et alJl2005^ . 
Table 1 lists the HOD parameters for the five values of erg 
considered in Paper I, with all masses scaled to Qm = 0.3. 
These parameters were chosen such that ^R{r) was nearly 
identical for all values of erg. We make mock galaxy distribu- 
tions by populating the N-body halos with galaxies accord- 
ing to equation 13511 . Satellite galaxies are placed randomly 
throughout the halo following the appropriate NFW profile 
for each mass and c osmology, calculated with the model of 
iBullock et alJ i200j) . Velocities are selected randomly from 
a Gaussian distribution with dispersion given by equation 
0. See Paper I for further details. 

Figure |H| plots the three redshift-space diagnostics for 
the central model. In all panels, the filled circles represent 
the N-body data, the solid line plots the model calculation, 
and the dotted line represents the dispersion model (eq. [5|)3 
for /3 = 0.46 and the best-fit value of a^ of 418 kms~i. Note 



2 Af jnin is set by the space density of galaxies once Mi and Osat 
are chosen. 

3 We use the implement ation of the dispersion model in configu- 
ration space as shown in lHawkins et alj i2003r . 
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20 
r [h"' Mpc] 

Figure 8. One-dimensional measures of redsliift-space cluster- 
ing. In each panel, the filled circles show the results of HODS 
from the N-body simulations (see Table 1). Error bars are shown, 
but at most r are smaller than the point size. Solid lines show the 
model calculations, and dotted lines show the best-fit dispersion 
model with dv = 418 kms~^. Panel (a) plots r^/2 against r^. 
Panel (b) plots fo/H against r. Panel (c) plots Qj against r. 



that our model does not have a parameter (like o-fc) adjusted 
to match the redshift-space distortions, and in that sense 
has less freedom to reproduce the data, though the formulas 
describing the halo pairwise velocities were calibrated on the 
simulations. Error bars on the N-body results are calculated 
from the dispersion of the five simulations divided by \/4 to 
yield the error in the mean, but they are generally too small 
to be seen outside the plot symbols. 

FigurelHK plots the results for r^/2 a-s a function of trans- 
verse separation. Both the N-body results and the model 
show a wave pattern that results from the transition in (,R{r) 
from the one-halo to two-halo term (see §5 in Paper I). As 
the dominant source of pairs shifts from one-halo pairs in 
high-mass systems to two-halo pairs in low-mass halos, the 
galax y velocity dispersion decreases sharply (Slosar ct alj 
[20061: iTinker et alj|2006i3) . The dispersion model shows no 
such wave pattern and is a poor fit to the N-body data, 
demonstrating the inability of a simple, scale-independent 
velocity parameter to describe small-scale galaxy motions. 

Figure IHJj plots ^o/h ^s a function of log(r). The large- 
scale infall amplifies the clustering, making ^o/r > 1 at r > 4 
h~^Mpc, while non-linearities suppress clustering, driving 
Co/Cfl to zero at small separation. Because Co('") is the low- 
est order multipole, it is least sensitive to the higher order 



moments of the velocity field. Thus, both the HOD model 
and the dispersion model provide good fits to the N-body 
results. 

Figure |HJ; plots Q^ as a function of r. Non-linearities 
become apparent in the quadrupole at a larger scale than in 
the monopole; although Q^ fiattens to a horizontal asymp- 
tote at large scales, the data deviate significantly from this 
asymptote at r < 20 h~^Mpc, and Qc becomes negative 
at r ~ 6 ^~^Mpc. Both the HOD model and the dispersion 
model reproduce the large-scale behavior, but the dispersion 
solution is a poor fit at r < 10 /i~^Mpc. The HOD model is 
accurate at all scales. 

Figure |^ plots the results of the N-body simulations 
for these same diagnostics but with variations in the 
{Q.m, erg, Qu) parameter space. In each panel, the circles plot 
the N-body data, with filled and open circles representing 
the high and low values of the isolated parameter, respec- 
tively. The solid and dotted lines are the model results for 
the high and low parameter values, respectively. The bot- 
tom row shows the results for the highest and lowest values 
of (Tg considered in Paper I (0.6 and 0.95), with Qm ~ 0.3 
and Ov = 1.0. The middle row presents the extreme values 
of Qm (0.1 and 0.5), with erg — 0.8 and Qf„ — 1. The top row 
presents results for setting Ov to 0.8 and 1.2 while keeping 
rim and (Tg fixed at their central values. In all cases, the 
HOD model accurately tracks the changes in the redshift- 
space diagnostic with changes in cosmological parameters, 
despite having no adjustable parameters of its own. 



4.2 Parameter Recovery 

In Paper I we concluded that combining redshift-space in- 
formation from multiple scales can break degeneracies and 
allow determination of flm, o"g, and «„. Models along the 
linear theory degeneracy axis, Q'^^'ag = const., have nearly 
identical large-scale redshift distortions that can mask sig- 
nificant changes in cosmology. Models along the degeneracy 
axis controlling satellite galaxy velocities, fimaS = const., 
have the same value of r^/2 when measured at r^ < 0.2 
h~^Mpc. Thus it is possible for widely varying cosmologies 
to produce the same large-scale anisotropies and the same 
small-scale FOGs. Despite these degeneracies, information 
from intermediate scales, where the one-halo to two-halo 
transition occurs, is able to break these degeneracies. 

To test the ability of the model to recover the correct 
cosmological parameters, we use this set of large scale/small 
scale degenerate models as the input data. The values of 
Jim and av for each model are listed in columns 6 and 7 of 
Table 1. The models are scaled to values of flm such that 
fH — 0.46 for all five HODs, and values of a„ are chosen such 
that rj/2 ~ 2.1 h~^Mpc at r^ = 0.1 h~^Mpc. All models 
have approximately the same real-space correlation function 
(see Figure 2 in Paper I). We use four observables; ^air), 
Co/fl) Qsi fi'iid r^/2, with diagonal error bars calculated from 
the standard deviation of the mean of the five simulations. 
We also add a systematic error to the model of 3% in each 
quantity to prevent data points with small errors from dom- 
inating the x^ minimization. With this systematic error, the 
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Figure 9. Measures of redshift-space clustering for various cosmologies. In each panel, the circles represent the N-body results and the 
lines show the model results. In each row, a cosmological parameter is varied between the high and low values used in Paper I. The 
top row plots results for a„ = 0.8 (open circles; dotted lines), and a^ = 1.2 (filled circles; solid lines). The middle row plots results for 
f2m = 0.1 (open circles; dotted lines), and Hm = 0.5 (filled circles; solid lines). The bottom row plots results for erg = 0.6 (open circles; 
dotted lines), and erg = 0.95 (filled circles; solid lines). 



X per degree of freedom of the model calculation with the 
N-body data is approximately one. 

Figures HUt i and 1 10b plot rj/2 and Q^, respectively, for 
the input N-body data. In panel (a), the models produce 
degenerate r^j^ values inside the one-halo term, Va < 0.7 
h~^Mpc. If av were set to unity for all five models, the 
r^/2 curves would separate at small radii due to the higher 
velocity scale of satellite galaxies for the higher Qm models. 
Outside the one-halo term, the models separate modestly, 
but are distinguishable outside their errors. In panel (b), all 
models have the same asymptotic value of Qj , but all models 
are below the linear theory prediction of 0.55 at all r for 
/3 = 0.46. The different values of a„ for each model cause 
the radii at which non-linearities set in to differ; for models 
with lower values of a„, Q( becomes negative at smaller r. 

To explore the likelihood of cosmological parameters 
from these data we us e the Monte Carlo M arkov chain 
(MCMC) technique (e.g.. lDoran fc Mulleil2o'o3) . In MCMC, 



a randomly chosen point in parameter space is automati- 
cally accepted as part of the chain if the x^ for that point 
is less than that of the previous point. If the x^ is larger, 
the probability of acceptance is exp(— Ax^/2). To choose the 
next point in parameter space, parameters are chosen from 
a Gaussian distribution centered on the previously accepted 
point. The widths of the Gaussians are the eigenvalues of 
the covariance matrix of the points already in the chain. 
Each chain is che cked for conver g ence u sing the power spec- 
trum method of iDunklev et alJ (|2005|). The benefit of the 
MCMC approach is that it allows one to marginalize over 
the HOD parameters; because different HODs will produce 
different galaxy velocity fields, it is necessary to account 
for uncertainties in the HOD parameters from the errors in 
the real-space correlation function. For each element in the 
chain, the total x^ is the sum of the x^ values for each of 
the four observables. 

For this analysis we have altered the parameters of the 
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Figure 10. Recovery of input cosmological parameters from MCMC analysis for the large-scale/small-scale degenerate HOD models 
with the parameter values listed in Table 1. Panels (a) and (b) plot the input N-body data for two of the four observables used. Panel 
(a) plots r^ /2 as a function of r^ for all five models. This quantity demonstrates the small-scale degeneracy of the models. Panel (b) 
plots Q^ as a function of r. This quantity demonstrates the large-scale degeneracy of the models. Panel (c) plots the 1- and 2-cr error 
contours in the Qm-crg plane, marginalizing over a„ and HOD parameters. The filled circles represent the input parameters. Panel (d) 
shows the error contours in the Qm-ctv plane, marginalizing over erg and HOD parameters. 



Jenkins mass function to better match the mass function 
of our own simulations. Due to shot noise in the particle 
distribution, the friends-of-friends algorithm systematically 
overestimates the masses of halo s with A'^ less than a few 
hundred particles ijWarren et al.l2 fl05'). At fixed (true) mass, 
this leads to an overestimate of th e mass function by a few 
percent. As shown in Appendix B oflTinker et al.! 1200 ST . the 
slight difference between our mass function and the Jenkins 
fitting function creates a discrepancy between the numerical 
and analytic one-halo terms of ^R{r). In the MCMC analy- 
sis, this discrepancy can drive the chain to the wrong solu- 
tion even though the difference in the model results is small. 



We stress that this modification is not a consequence of the 
model; the model will correctly describe the redshift-space 
distortions for any input mass function. The parameters of 
the Jenkins formula are (0.315,0.61,3.8). Our modified pa- 
rameters are (0.345, 0.636, 3.43), which increases the number 
of halos less than M, by ~ 7%. The use of simulations with 
higher mass resolution would prod uce mass funct i ons in bet- 
ter agreement with the results of I Jenkins et al.l l|2001il and 
ImrreniraOiooi). 

Figures IT?7t and llOH show the constraints on the cosmo- 
logical parameters from the MCMC results. Panel (c) shows 
the ilm-crg projection of the parameter space, and panel (b) 
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Figure 11. Constraints on f2m, Tgi and Ot, for models with extra free parameters. In each panel, the dotted contours represent the error 
contours from Figure Hni for HOD3. The top three panels show constraints in the ilm-oiv plane, while the bottom panels show constraints 
in the n^-crg plane. In panels (a) and (b), the solid lines correspond to a model in which F is introduced as a free parameter, in addition 
to the HOD and cosmological parameters. In panels (c) and (d), the solid lines correspond to a model in which the ratio between the 
halo concentrations and the concentration parameter for satellite galaxies, Cg^i/cj^g^iQ is a free parameter. In panels (e) and (f), the solid 
lines correspond to a model in which {Nsa.t)M contains an extra parameter Mcut (see equation I36I V 



panel rotates the data to show the Q,m-av projection. The 
solid circles plot the input cosmological values. The con- 
tours show the one- and two-a confidence levels on each 
parameter. For most of the models listed in Table 1, the 
correct values are within either the one- or two-cr error con- 
tour. For the central model, the MCMC analysis correctly 
yields O^ = 0.306 ± 0.006 and ag = 0.799 ± 0.009, but the 
velocity bias q„ = 0.955 ± 0.015 is 3-(t ofT the true value 
of 1.0. For models HODl, HOD2, and H0D4, the analytic 
model correctly recovers the true a^ to within 1-a (~ 1.5%). 
The source of this discrepency in HODS can be traced to 
^H(r); the analytic model underpredicts the amplitude of 
^ fl(r) for this model at r « 1 /i~^Mpc (see Appendix B 
in lTinker et al.ll2005l for a comparison between the analytic 
model and N-body results for this HOD). To produce a bet- 
ter fit to the real-space clustering, a higher value of Qsat is fa- 
vored (1.032 versus 1.005). This increases the mean velocity 
dispersion of galaxies, so the likelihood is higher for a model 
with slightly negative velocity bias, ~ 4.5% below the true 
value. This discrepancy is not significant in model the other 
models (we will discuss the results for HOD5 presently). In 
tests with N-body HOD models tailored toward samples of 



brighter, less abundant galaxies, this discrepancy is absent 
in the analytic model for £,R{r). 

Compared to HOD models 1 through 4, which recover 
the input cosmology to better than ~ 2% (with the excep- 
tion discussed in the previous paragraph), the recovered pa- 
rameters for H0D5 are barely within the 2-a contours. The 
analysis yields Q.m = 0.491 ± 0.014, erg = 0.607 ± 0.010, 
and a„ = 0.765 ± 0.012, compared with the input values of 
(0.47, 0.6, 0.797). The high preferred value of fi™ drives the 
preferred value of Qt, lower than the input value. The most 
likely value of ag, however, is unaffected. The larger errors 
for this model could represent a breakdown in O!^ as an 
approximation to the linear theory scaling of halo velocities 
or the inability of linear theory itself to accurately model 
the scaling. Further testing is required to isolate the source 
of the error, although it is likely outside the range of interest 
for the true value of Q.m ■ 



4.3 Adding freedom to the model 

The test in the previous section made several simplifying 
assumptions. First, the shape of Piin(fc) was fixed to the 
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true value from the simulations. Second, the radial profiles of 
the satellite galaxies were assumed to follow the true profiles 
used to populate the simulations. Third, the functional form 
of the HOD used in the analysis was the same as that used 
to create the test data. 

In Figure 1111 we test the robustness of the analytic 
model against these assumptions. Panels (a) and (b) com- 
pare the results of the MCMC analysis from Figure 1101 to 
an analysis in which the power spectrum shape parameter F 
is added as an additional free parameter with no prior. The 
dotted contours represent the 1- and 2-a constraints on the 
cosmological parameters from Figure 1101 in which F is set 
to 0.2. The solid contours are the results once marginalized 
over F. The constraints are nearly identical, even with the 
added freedom to the model, and the most likely value of 
F = 0.204 ± 0.005. As discussed in §1, altering the shape 
of -Piin(fc) has limited impact on the redshift-space observ- 
ables employed here. But altering F changes the amplitude 
of S,R{r) for a fixed HOD, so large deviations from the true 
value of F are easily excluded by the real-space clustering 
alone. 

Figures UTt and lllH show the cosmological constraints 
of the analytic model when the galaxy density profile within 
halos is allowed to vary. In this model, we introduce a new 
parameter a^ = Cgai/chaio, where c represents the concen- 
tration par ameter of the NFW density profile, and Chaio is 
taken from iBuUock et alj 1120011) . This parameter preserves 
the shape of the dependence of c on halo mass, but allows its 
normalization to vary. Panel (c) presents the constraints in 
the Qm-Qv plane from Figure HUl and from varying-ttc model. 
Marginalizing over a^ has a noticeable efi'ect on the con- 
straints. The 1-a errors on f2m increase from 0.006 to 0.009, 
while the 1-a errors on a„ increase from 0.015 to 0.021. Al- 
though the most likely values of ilm and a^ are consistent 
in the solid and dotted contours, marginalizing over Oc does 
not expand the contours in a symmetric fashion. As dis- 
cussed in the previous section, the analytic model for HODS 
underpredicts £,R{r) at r = 1 h~^Mpc. When marginalizing 
over Uc, this discrepancy with £,R{r) can be compensated 
for by decreasing the concentrations of the satellite galax- 
ies, thereby shifting one-halo pairs out to larger separation. 
Thus, in the varying-Oc analysis, the preferred value of Qsat 
is closer to the true value of 1.005, skewing the error con- 
tours, but the constraints on Osat are somewhat weaker than 
in the fixed-Oc analysis, yielding larger overall errors. The 1- 
(j constraints on Ac are 0.93 ± 0.16. Panel (d) presents the 
constraints in the flm-o's plane for these two models. The 
1-0" error contours are the same in the two models, but the 
extra freedom of the varying-flc model produces a small ex- 
tension of the 2-a contour down the Q'^ag degeneracy axis. 
The weaker constraints produced by the varying-Oc model 
are apparent in Figure ITTl but they only increase the errors 
from ~ 2% to ~ 3% on the cosmological parameters. 

Lastly, we demonstrate the robustness of the model to 
extra freedom in the parameterization of the HOD itself. 
The satellite occupation function used throughout this pa- 
per is a simple power law with two free parameters. In Fig- 
ures lllb and lilt we compare the MCMC results from our 
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Figure 12. Left-hand panels plot the standard redshift-space 
observables for HOD parameters constrained to match the Mr — 
Slog h < —21 SDSS Wp(rp). The filled circles are N-body results 
for r^/2 in panel (a), ^o/H ™ panel (b), and Qj in panel (c). The 
solid lines are the analytic model calculations for each quantity 
for the same HOD parameters. Right-hand panels plot the same 
observables, but for a model in which F = 0.12. The cosmology 
and HOD parameters are the same as HODS in Table 1. Points 
and lines are the same as panel (a)— (c). 



two-parameter {Nsat)M to a model with a three-parameter 
{A'"sat)A/ of the form 



{iVsat)A/ = 



M - Mc, 
Ml 



(36) 



where Mcut is a cutoff mass scale. In equation 13511 . {A'sat)Af 
is cut off sharply at Mmin- The form of equation I36II has a 
soft cutoff that can be at any mass, although halos below 
Afmin are not allowed to host satellites. (Even with this re- 
striction, models with Mcut < Mmin are not degenerate until 
Mcut ^ Mmin.) However, this extra freedom in the occupa- 
tion function produces a negligible effect on the cosmological 
constraints. Panels (e) and (f) once again compare the error 
contours in Qm-ctv and Qm-as planes, respectively, for the 
standard two-parameter (A'^sat)*/ and varying-Mcut model. 
The contours are nearly unchanged between the two models. 



4.4 Other mock samples 

The HOD parameters in Table 1 produce galaxy space den- 
sitie s consistent with t he Mr — 5 log h < —20 SDSS sam- 
ple JZehavi et alJl2005h . Although our N-body simulations 
lack the resolution to model fainter samples, we can test 
the accuracy of the model against mock data matching to 
brighter samples. In Figures ll2M12b . we compare the an- 
alytic model to N-body data created with HOD parame- 
ters determined from fitting Wr,(r „) measurem e nts o f the 
Mr - 5 log ft < -21 sample from IZehavi et al.l ll2005l) . To 
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populate the si mulations, we use th e HOD parameters listed 
in Table 3 of IZehavi et^ (|200,4 l: Mmi„ = 5.25 x 10^^ 
ft-^M©, Ml = 12.3 X 10" ft-^Mo, and a^at = 1.39. These 
HOD parameters are used to populate the simulations with 
as = 0.9, and the halo velocities are scaled to fim = 0.3. 
This sample yields a space density lower than the HODs 
in Table 1 by a factor of ~ 5, and emphasizes higher halo 
masses. The N-body results are shown with the filled cir- 
cles for rj/2 in panel (a), for £,o/r in panel (b), and Q^ in 
panel (c). The model calculations for the same cosmology 
and HOD are shown with the solid lines in each panel. For 
r^/2, the lower space density of galaxies increases the error 
on this diagnostic within the one-halo term, partially from 
Poisson noise and partially from a less robust determination 
of the asymptotic value of ^(r^) at small r-^. But the ana- 
lytic model compares favorably to the data, correctly pre- 
dicting the location and depth of the minimum at r,,- = 1.5 
h~^Mpc. The analytic model also accurately predicts the 
shape of ^q/r, which has a sharper transition from linear to 
non-linear scales than the results in Figure |HJ) due to the 
emphasis on higher mass halos. For Q^, the analytic model 
correctly describes the large-scale asymptote as well as the 
scale at which non-linearities set in. 

In Figures [T^ - I12L we compare the analytic model to 
N-body data created with the HOD parameters and cosmol- 
ogy of the central model, but populating the halos of sim- 
ulations in which the shape parameter of the linear matter 
power spectrum is F = 0.12.* Although a power spectrum of 
this shape is unable to yield acceptable HOD fits to current 
SDSS measurements of Wp{rp) (Z. Zheng, private commu- 
nication), it is worth testing the analytic model against an 
extreme change in Piin(fc). Changing the shape of Piin(fc) al- 
ters not only the halo mass function and bias relation but 
the one-point distribution of dark matter densities used in 
the velocity model as well. For all three observables, the an- 
alytic model performs well in modeling the redshift-space 
clustering for a substantially different power spectrum. 

This success is a strong indication that the analytic 
model is correctly describing the underlying physics of 
redshift-space distortions as it accurately describes a cos- 
mological model quite different from the one on which it 
was calibrated. 



5 SUMMARY 

We have presented a model for the redshift-space correlation 
function of galaxies using the framework of the halo occupa- 
tion distribution (HOD). The HOD quantifies galaxy two- 
point clustering by separating the total number of galaxy 
pairs into the contributions from within a single halo and 
pairs from two distinct halos. In redshift space, the one-halo 
term is obtained by assuming galaxies follow an isother- 
mal, Gaussian velocity distribution with a dispersion pro- 
portional to the virial dispersion of the dark matter halo. 



* The specifics of these simulations are identical to the F 
simulations. See Paper I for more discussion. 



0.2 



The constant of proportionality Qf„ represents the veloc- 
ity bias of galaxies, and it encompasses the effects of non- 
isothermality, anisotropy, or systems being unrelaxed. Cen- 
tral galaxies are assumed to follow the center of mass of the 
halo they occupy. Pairs of galaxies from two distinct halos 
move with a combination of virial motions and the halo cen- 
ter of mass velocity. Halo velocity statistics have significant 
third and fourth moments which must be modeled to cal- 
culate redshift-space clustering. The key simplification that 
makes our analytic model tractable is that the pairwise dis- 
tribution of halo velocities at fixed density is approximately 
Gaussian and approximately independent of halo mass. The 
non-Gaussian shape of the full distribution arises from aver- 
aging over pairs in difi^erent density environments, and the 
dependence of the distribution on halo mass arises because 
higher mass halos tend to reside in denser environments. 

Our approach solves several complications in modeling 
pairwise velocity statistics at non-linear and quasi-linear 
scales. It demonstrates why halo velocity PDFs in simu- 
lations are distinctly non-Gaussian, with significant skew- 
ness and kurtosis that fluctuate with scale. The skewness 
arises from halo pairs in high-density en vironments, where 
pairs have high relative velocities (see also iJuszkiewicz et alJ 
1998!). The kurtosis results from a superposition of Gaus- 
sians of varying widths from halo pairs in wide-ranging en- 
vironments. The model also circumvents to problem of cor- 
related radial and tangential velocities. At a given 5, Vr and 
vt are independent, and it is straightforward to calculate Vz 
at any angle with respect to the observer's line of sight. The 
model is highly flexible in that it can describe redshift-space 
clustering for samples of any space density without the asso- 
ciated complications arising from N-body simulations that 
sacriflce volume for mass resolution, or resolution in favor of 
box size. It also eliminates the need to run large numbers of 
simulations when exploring parameter space. 

This model can be used to obtain constraints in 
the {Q,rn,<ys,av) parameter space from redshift-space 
anisotropies following the blueprint detailed in Paper I. Ob- 
servations of Wp{rp) provide tight constraints on the pa- 
rameters of the HOD, separating the constraints on bias 
from the constraints on cosmology. Our choices for redshift- 
space observables largely normalize out the shapes of the 
real-space galaxy correlation function and the linear matter 
power spectrum. Marginalization over Piin(fc) and HOD pa- 
rameters do not degrade the ability of the model to constrain 
parameters. The only unknown quantity that notably affects 
these constraints is the proflle of satellite galaxies within 
groups Cgai, a quantity that can be observationally mea- 
sured ([Y ang et all l2005fl and theoretically pr edicted (e.g., 
iDiemand et alj|2004 iNagai fc Kravtsovl |2005^ . Even with 
no prior on Cgai, cosmological constraints only increase by 
~ 1% relative to the true value. 

Our HOD approach represents a significant improve- 
ment over the dispersion model (Eq. |2]) for two main rea- 
sons. First, the HOD model circumvents the use of the 
redshift-space distortion parameter /3 entirely. The cosmo- 
logical parameters Qm and as are direct inputs of the model, 
eliminating the degeneracy involved with using /3 to deter- 
mine the matter density. Second, there are no free parame- 
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ters like at, which creates a degeneracy with (3 as well. The 
velocity bias parameter Qt, is a truly physical parameter that 
is cosmologically relevant and is a testable prediction of nu- 
merical simulations of galaxy formation and semi-analytic 
models. Measurements of q:„ can reduce uncertainty in dy- 
namical mass estimates of galaxy groups, utilized in group 
catalogs culled f rom large-s cale optical surveys (e.g., the 
2PIGG catalog of lEke et elk2004.) . Alternately, studies that 
use the abundance of systems as a function of their velocity 
dispersion require knowledge of ap to provide useful compar- 
isons to theoretical models f,Dayis et al. 2005; Gcrkc ct al. 
120051) . In this paper we have assumed that central galaxies 
are not subject to virial motions. Modest values of a^c have 
limited effect on redshift-space clustering, but non-zero a^c 
can easily be incorporated into the analysis of observational 
data and either marginalized over or constrained if the data 
are of sufBcient precision. 

Utilizing clustering information on multiple scales 
breaks the two primary degeneracy axes in redshift-space 
distortions. At large scales, linear theory correctly predicts 
that models with the same value of /3 will have the same 
large-scale distortions (even though linear theory is a poor 
descriptor of the shape of the distortions for a given value 
of (3). This creates the degeneracy Q^ag = const. At 
small scales, non-linearities are determined by the veloc- 
ity scale of dark matter halos, creating the degeneracy axis 
ClrnOp = const. Paper I concluded that intermediate-scale 
clustering, 1 /i~^Mpc < r < 10 /i~^Mpc, provides the key 
for breaking the large-scale/small-scale degeneracies. Here 
we have quantitatively demonstrated this result. The typ- 
ical "banana-curve" degeneracy axes can be seen in Fig- 
ures nut and llOt l. but the simulations are clearly distin- 
guished from one another by the analytic model. The large- 
scale distortions set the amplitude of the ^^as degeneracy 
curve, the small-scale distortions normalize the f2ma^ de- 
generacy curve, and the intermediate-scale clustering locates 
the model along both curves. 

The current data releases of the SDSS and 2dFGRS 
have reached the necessary volume and statistics to imple- 
ment this HOD analysis. The ability to extract information 
from galaxy clustering measurements at quasi-linear and 
non-linear scales enables precise constraints on our speci- 
fied subset of cosmological parameters that can be compet- 
itive with large-scale measures such as CMB anisotropics, 
Lyman-Q forest, and the large-scale galaxy power spectrum. 
Although constructed on more fundamental physics, these 
methods are sensitive to quantities like the CMB tensor-to- 
scalar ratio, mass of the neutrino, and curvature of the in- 
flationary fluctuation spectrum. For small-scale galaxy clus- 
tering, modest changes in the shape of Piin(fc) have marginal 
impact on halo occupation constraints. Thus redshift-space 
clustering measurements provide a powerful and comple- 
mentary method for estimating the amplitude of dark mat- 
ter clustering and the mean matter density of the universe. 
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APPENDIX A: 
TERM 



OUTLINE OF TWO-HALO 



For the convenience of readers who wish to implement our 
model, we present a concise outline of the calculation of 
£,2h{rc,rTy), beginning with the definition of the streaming 
model. 



l + i2Y,(r„ 



[l + i2^{r)]P2^(vAr,4>)dv. 



(Al) 



where r^ — r^ + z^ , v^ — H {r^r — z), and cos(f) — Tajr. 
C2h('') is the two-halo contribution to the real-space corre- 
lation function, and P2h(t'z|'', 0) is the line-of-sight velocity 
PDF of galaxy pairs from two distinct halos. The velocity 
PDF in equation lAH is a pair-weighted average over all 
possible combinations of halo masses. For a given halo pair, 
M\ and M2, the line-of-sight PDF is an integral over all 
possible densities in which that pair can exist. 

The line-of-sight PDF of a halo pair at density & is 



P9+h(t;^|<^,r,Mi,M2,5) = 



4 

E 

1=1 



l-KOi 



■ exp 



-{Vz 



fir sm ( 



2a? 



(A2) 



where each term in the summation represents a different 
combination of central and satellite galaxy pairs from the 
two halos. The dispersion for each term is 



/ 2 2 
(a* cos I 



2 . 2 
+ Ur sm ( 



+ k 



(i)i2 
gal J ' 



(A3) 



where at and Ur are the tangential and radial dispersions 
of the dark matter halos, and a^J^^ is the dispersion of the 
galaxy pair from internal halo motions, defined by equation 
0. The radial dispersion is 



ar = 200f" 
where 



0.3/ V0.8/ \p2m.r 



(A4) 



p200,r(Ml,M2,r) 



5.0i?i/,'i 



+ 



11. 5R 



1/2 



and 

a(r) 



35/i-iMpc 



+ 0.50, (A5) 



(A6) 



© 2005 RAS, MNRAS 000, 000-000 



20 J. L. Tinker 



For the tangential velocities, p200,r is replaced in equation 
ljA4|l with 



P200,t{Mi,M2,r) = 



Miim,i('-) = ■r-_Rvir(Mmin) 
Mii„,,2(r) = r-_Rvir(Afi). 



(A14) 



Equation 1A13I I is inserted directly into equation lAll . 



,1/2 



+ 



,1/2 



7.2i?X:i / V 12.6i?tir,o , 



+ 0.48. (A7) 



The units of r and i?vir are h ^Mpc. The mean velocity 
fj, of equation IIA2l l is given by the spherical collapse model 
(eq. [IHl) at r ^ 4 h^^Mpc, linear theory (eq. HKl'l aX r > 2Q 
/i^^Mpc, and a linear combination of the two at intermediate 
scales, 



^(r, 8) =W ^sc + (1 — Ul) /Xlin, 

with 

( 1 if r :^ 4 

w= \ 1.86 - 0.62 In r if 4 sC r < 20 
I if r < 20 . 



(A8) 



(A9) 



The distribution function Pgj^h{vz\r, (j), Mi, M2, 5) is in- 
tegrated over all 5 to obtain 

Pg+h{v,\(t),r,Mi,M2) = 

Pg+hiv4(l,,r,Mi,M2,S)Pr„{S\r,Mi,M2)dS, (AlO) 



where the PDF of dark matter density used in equation 
JAlOl l is given by 



Pn{S\r,Mi,M2) = 

"l.45(&i+b2) + (r/9.4i?vir,i)'''-' 



j4exp 



and 



Pm{S\r) 



2Tra', 



exp 



1 + 5 



[ln(l + S) + G\l2f 



Pm[&\r) 



2o{ 



(All) 



(A12) 



where cri(r) = ln[l + (T^(r)] and am(r) is the mass variance 
in top-hat spheres of radius r. To calc ulate Omir^, we use 
the non-linear matter power spectrum of lSmith et al.l (120031) 
to better describe the distribution of dark matter at small 
scales. To calculat e the halo bias factors, we use the function 
in Appendix A of lTinker et alJ (|2005|). The final step is to 
integrate over all halo pairs, 

/•Mlim 1 J 

P2h(^;.|r,0) = 2(4)-2 / ' dMi-^{iV)Mi 



■ dMi 

"""'" dn 

rfM2^^{7V)M,Pg+h(t'.|r,0,Mi,M2) , (A13) 



where fig is given by equation @ , and the limits of the inte- 
gral are determined by the condition that no halo pair can 
be closer than the sum of their virial radii. Using the rela- 
tion 7?vir = (3M/47rpAvir)^ (recall that we have defined 
Avir = 200 throughout this paper). 



REFERENCES 



Abazajian, K., et al. ApJ (submitted), |astrc>-p h/0408003] 
Bean, A. J., Ellis, R. S., Shanks, T., Efstathiou, G., & 

Peterson, B. A. 1983, MNRAS, 205, 605 
Benson, A. J. 2001, MNRAS, 325, 1039 
Berlind, A. A., & Weinberg, D. H. 2002, ApJ, 575, 587 
Berlind, A. A., Weinberg, D. H., Benson, A. J., Baugh, C. 

M., Cole, S., Dave, R., Frenk, C. S., Katz, A., & Lacey, 

C. G. 2003, ApJ, 593, 1 
Bond, J. R., Cole, S., Efstathiou, G., & Kaiser, N. 1991, 

ApJ, 379, 440 
Bullock, J. S., Kolatt, T. S., Sigad, Y., SomerviUe, R. S., 

Klypin, A. A., Primack, J. R., Dekel, A. 2001, MNRAS, 

321, 559 
Cole, S., Fisher, K. B., & Weinberg, D. H. 1994, MNRAS, 

267, 785 
Cole, S., Fisher, K. B., & Weinberg, D. H. 1995, MNRAS, 

275, 575 
Coles, P., & Jones, B. 1991, MNRAS, 248, 1 
CoUess, M. et. al. 2001, MNRAS, 328, 1039 
Collister, A. A., & Lahav, O. 2004, MNRAS (submitted), 

astro-ph/0412516 
Cooray, A. & Sheth, R. K. 2002, Phys. Rep., 372, 1 
Cooray, A., 2004, MNRAS, 348, 250 
Davis, M. & Peebles, P. J. E. 1983, ApJ, 267, 465 
Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 

1985, ApJ, 292, 371 
Davis, M., Gerke, B. F., & Newman, J. A. 2005, ASP Con - 

ference 339: Observing Dark Energy, fastro-ph/0408344) 
Diemand, J., Moore, B., & Stadel, J. 2004, MNRAS, 352, 

535 
Doran, M. & Miiller C. M. 2004, JCAP, 409, 3 
Dunkley, J., Bucher, M., Ferreira, P. G., Moodley, K., & 

Skordis, C. 2005, MNRAS, 356, 925 
Efstathiou, G., Bond, J. R., & White, S. D. M. 1992, MN- 
RAS, 258, 1 
Eke, V. R., et al. 2004, MNRAS, 348, 866 
Faltenbacher, A., Kravtsov, A. V., Nagai, D., & Gottlober, 

S. 2004, MNRAS (submitted), astro-ph/0408488 
Fisher, K. B., Davis, M., Strauss, M. A., Yahil, A., & 

Huchra, J. P. 1994, MNRAS, 267, 927 
Fisher, K. B. 1995, ApJ, 448, 494 
Fisher, K. B. & Nusser, A 1996, MNRAS, 279, LI 
Gerke, B. F. et al. 2005, ApJ, 625, 6 
Ghigna, S., Moore, B., Governato, F., Lake, G., Quinn, T., 

& Stadel J. 2000, ApJ, 544, 616 
Hamana, T., Issha, K., Yoshida, N., Suto, Y., & Jing., Y. 

P. 2003 MNRAS, 343 1312 
Hatton, S. & Cole, S. 1999, MNRAS, 310, 1137 
Hawkins, E., Maddox, S., Cole, S., et al. 2003, MNRAS, 

346, 78 
Jenkins, A., Frenk, C. S., White, S. D. M., Colberg, J. M., 



© 2005 RAS, MNRAS 000, 000-000 



Cole, S., Evrard, A. E., Couchman, H. M. P., & Yoshida, 

N. 2001, MNRAS, 321, 372 
Jing, Y. P., Mo, H. J., & Boerner, G. 1998, ApJ, 494, 1 
Juszkiewicz, R., Fisher, K. B., Szapudi, I. 1998, Apj, 504, 

LI 
Juszkiewicz, R., Springel, V., Durrer, R. 1999, Apj, 518, 

L25 
Kaiser, N. 1987, MNRAS, 227, 1 

Kang, X., Jing, Y. P., Mo, H. J., & Borner, G. 2002 MN- 
RAS, 336, 892 
Klypin, A., Gottlober, S., Kravtsov, A. V., and Kohklov, 

A. M. 1999, ApJ, 516, 530 
Kofman, L., Bertschinger, E, Gelb, J. M., Nusser, A., & 

Dekel, A. 1994, Apj, 420, 44 
Kravtsov, A. V., Berlind, A. A., Wechsler, R. H., Klypin, 

A. A., Gottlober, S., AUgood, B., Primack, J. R. 2004, 

Apj, 609, 35 
Lin, Y.-T., Mohr, J. J., & Stanford, S. A. 2004, ApJ, 610, 

745 
Mo, H. J., & White, S. D. M. 1996, MNRAS, 282, 1096 
Ma, C. & Fry, J. N. 2000, ApJ, 543, 503 
Madgwick, D. S., et al. 2003, MNRAS, 344, 874 
Nagai, D. & Kravtsov, A. V. 2005, ApJ, 618, 557 
Park, C, Vogeley, M. S., Geller, M. J., & Huchra, J. P. 

1994, ApJ, 431, 569 
Pavlidou, V. & Fields, B. D. 2005, PhysRevD, 71, 043510 
Peacock, J. A. & Dodds, S. 1994, MNRAS, 267, 1020 
Peacock, J. A., Smith, R. E. 2000, MNRAS, 318, 1144 
Peacock, J. 1999, Cosmological Physics, Cambridge Uni- 
versity Press 
Peebles, P. J. E. 1976, ApJ, 205, L109 
Peebles, P. J. E. 1979, AJ, 84, 730 
Peebles, P. J. E. 1980, The Large-Scale Structure of the 

Universe 
Percival, W. J. et al. 2002, MNRAS, 337, 1068 
Sargent, W. L. W. & Turner, E. L. 1977, ApJ, 212 L3 
Scoccimarro, R., Sheth, R. K., Hui, L., & Jain, B. 2001, 

ApJ, 546, 20 
Scoccimarro, R. 2004, PhysRevD, 70, 063001 
Seljak, U. 2000, MNRAS, 318, 203 
Seljak, U. 2001, MNRAS, 325, 1359 
Sheth, R. K. 1996, MNRAS, 279, 1310 
Sheth, R. K. & Tormen, G. 1999, MNRAS, 308, 119 
Sheth, R. K. & Diaferio, A. 2001, MNRAS, 322, 901 
Sheth, R. K., Hui, L., Diaferio, A., & Scoccimarro, R. 

2001a, MNRAS, 325, 1288 
Sheth, R. K., Diaferio, A., Hui, L., & Scoccimarro, R. 

2001b, MNRAS, 326, 463 
Sheth, R. K. & Tormen, G. 2002, MNRAS, 329, 61 
Skibba, R., Sheth, R. K., Connolly, A. J., & Scranton, R. 

2005, MNRAS (submitted), astro-ph/0512463 
Slosar, A., Seljak, U., & Tasitsiomi, A. 2006, MNRAS, 366, 

1455 
Smith R. E., et al. 2003, MNRAS, 341, 1311 
Spergel, D. N., et. al. 2003, ApJ, 148, 175 
Springel, V., Yoshida, N., & White, S. D. M. 2001 NewA, 

6, 79 
Tegmark, M., et. al. 2004a, ApJ, 606, 702 
Tegmark, M., et. al. 2004b, PhysRevD, 69, 103501 



Redshift- Space Distortions with the HOD 21 

Tinker, J. L., Weinberg, D. H., Z heng, Z., & Zehavi, L 2005, 

ApJ, 631, 41 astro-ph/0411777 
Tinker, J. L., Weinberg, D. H.,& Zheng, Z., L 2006a, MN- 
RAS (accepted), astro-ph/0501029 
Tinker, J. L., Norberg, P., Weinberg, D. H., & Warren, 

M. S. 2006b, ApJ (submitted), astro-ph/0603543 
van den Bosch, F. C, Weinmann, S. M., Yang, X., Mo, 

H. J., Li, C, & Jing, Y. P., 2005, MNRAS (submitted), 

astro-ph/0502466 1 
Warren, M. S., Abazajian, K., Holz, D. E., & Teodoro, L. 

2005, ApJL, submitted, ( astro-p h/05063 95 ) 
White, M. 2001, MNRAS, 321, 1 
Yang, X., Mo, H. J., & van den Bosch, F. C. 2003, MNRAS, 

339, 1057 
Yang, X., Mo, H. J., Jing, Y. P., van den Bosch, F. C, & 

Chu, X. Q. 2004, MNRAS, 350, 1153 
Yang, X., Mo, H. J., van den Bosch, F. C, Weinmann, 

S. M., Li, C, & Jing, Y. P. 2005, MNRAS, submitted, 

(astro-ph/0504477) 
York, D., et al. 2000, AJ, 120, 1579 
Yoshikawa, K., Jing, Y. P., & Borner, G.2003, ApJ, 590, 

654 
Zehavi, L, Weinberg, D. H., Zheng, Z., Berhnd, A. A., Frie- 

man, J. A., et al. 2004, Apj, 608, 16 
Zehavi, I., et al. 2005, Apj, 630, 1 
Zheng, Z. 2004, Apj, 610, 61 
Zheng, Z., Tinker, J. L., Weinberg, D. H., & Berhnd, A. A. 

2002, Apj, 575, 617 
Zheng, Z., Berlind, A. A., Weinberg, D. H., Benson, A. J., 

Baugh, C. M., Cole, S., Dave, R., Frenk, C. S., Katz, A., 

& Lacey, C. G. 2005, Apj, 633, 791 
Zurek, W. H., Quinn, P. J., Salmon, J. K., & Warren, M. S. 

1994, ApJ, 431, 559 



© 2005 RAS, MNRAS 000, 000-000 



